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We present a formula to evaluate the spontaneous emission spectra of an atom in contact with a 
radiation field with non-Markovian effects. This formula is written in terms of a two-time correlation 
of system observables and the environmental correlation function, and depends on the distance 
between the emitting atom and the detector. As an example, we use it to analyze the fluorescence 
spectra of a two level atom placed as an impurity in a photonic crystal. The radiation field within 
those materials has a gap or discontinuity where electromagnetic modes cannot propagate in the 
stationary limit. In that situation, the atomic emission occurs in the form of evanescent waves which 
are detected with less efficiency the farther we place the detector. The methodology presented in 
this paper may be useful to study the non-Markovian dynamics of any quantum open system in 
linear interaction with a harmonic oscillator reservoir and within the weak coupling approximation. 



INTRODUCTION 



Most of the theoretical work developed in the context of quantum open systems relies in the use of the Markov 
approximation, vifhich assumes that when a quantum open system (i.e. the atom) interacts with its environment (for 
instance a radiation field), the later is not affected in a time scale comparable to that of the system. In other words, 
the radiation field absorbs the emitted photon instantaneously in comparison to the atomic evolution time scale. 

However, the Markovian assumption is only valid for radiation fields in which the photonic density of states (DOS) 
is a continuous and smoothly varying function in the frequency space. When this function varies in a range which 
is comparable to the atomic spontaneous emission rate F, then the environmental relaxation time (also correlation 
time) Tc, can no longer be considered negligible in comparison with the atomic relaxation time. We are then dealing 
with non-Markovian interactions, in which the environment does not absorb instantaneously every photon it receives 
from the atom, but needs a certain memory time to do it and return to its equilibrium state. 

Non-Markovian interactions may occur when the radiation field is placed in some structured environment, so that 
its photonic DOS may present sudden jumps or singular behavior in a frequency range comparable to the spontaneous 
emission rate. For instance, in the so-called photonic crystals (PC) [H, B], there is a periodicity in the refraction index 
which produces the scattering of light outside of the crystal for certain frequencies that are related to the lattice 
periodicity. Since these modes are absent from the crystal, the photonic DOS is zero in the corresponding frequencies, 
so that a gap is formed. The frequency ranges in which the photonic DOS is non-zero are known as bands. Within 
the bands, the photonic DOS varies in frequency range comparable to the atomic spontaneous emission rate, which 
gives rise to long memory time for the environment. This kind of non-Markovian interaction may also be found in 
a quantum cavity, when the range of variation of the photonic DOS is comparable to the spontaneous emission rate 
ili- 

In non-Markovian interactions the emitted photon may be eventually reabsorbed by the atom, which produces a 
back-action in the system dynamics, both in the evolution of its expectation values and the multiple-time correlation 
functions (MTCF) of its observables. An important consequence of this is that some related physical quantities such 
as the atomic emission spectra also display non-Markovian features. Hence, it is necessary to obtain a more general 
formula for the atomic emission spectra without using the Markovian approximation. 

In most of the literature (see for instance [3] and references therein), the analytical expression for the spectra is 
derived by assuming the Markov approximation at some stage. To be more specific, in photodetector experiments 
the emission spectra depends on a two-time correlation of the positive and negative operators of the radiation field 
emitted by the atom. In order to evaluate it, this two-time correlation is expressed in terms of a two-time correlation 
of system observables, which is done by assuming the Markovian approximation. In the end, it is find that within this 
approximation the spectra is related to a Fourier transform of the system two-time correlation 0, Moreover, once 
the two-time correlation is computed through the Quantum Regression Theorem, a Lorentzian profile for the spectra 
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is obtained. This theorem states that the evolution of A''-time correlation functions can be computed with the master 
equation by only considering a rotation in its initial condition. 

In this paper, we derive a formula to compute the atomic emission spectra for non-Markovian interactions. We 
show that the two-time correlation of the positive and negative operators of the radiation field is related to a two- 
time correlation of system observables through a double convolution integral, that also depends on the environmental 
correlation function. The decaying of this environmental correlation function gives the correlation time or relaxation 
time of the environment Tc, which for non-Markovian interactions is non-zero. Moreover, the two-time correlation of 
the system observables can no longer be described within the so-called Quantum Regression Theorem (QRT), since 
the derivation of such formula relies in the Markov approximation 0, S, 0|- Iii the non-Markovian case, the set of 
evolution equations derived in 13, 11 1 should be used to evaluate MTCF. Moreover, the formula we derive for the 
emission spectra is dependent on the spatial distance between the emitting atom and the detector. 

In order to illustrate the theory, we study the fluorescence spectra of a two-level atom in contact with the modified 
radiation field that exists within a photonic crystal. In the stationary limit, the emission in the frequency region 
corresponding to the gap is made in the form of evanescent modes, which present an exponentially decaying spatial 
dependency. It is clear that when evanescent modes are present in the emission, the fluorescence spectra depends 
dramatically on the distance from the atom to the detector. 

The plan of the paper is the following. In Sec. we present the general Hamiltonian for the class of systems that 
we want to describe, as well as the environmental correlation function that characterizes the interaction. In Section 
we present a formula to obtain a non-Markovian spectra in terms of a two-time correlation of system observables, 
convoluted with the environmental correlation function. This formula depends explicitly on the distance between the 
emitting system and the detector, and is one of the main results of this paper. Section is devoted to review the 
evolution equations of system two-time correlations within the weak coupling limit. In the last part of the paper 
(Section ) we apply the theory to compute the fluorescence spectra of a two-level atom placed as an impurity in a 
photonic crystal, showing that the distance between the atom and the detector is an important parameter to take 
into account if we want to characterize correctly the spectra in the gap region. 



NON-MARKOVIAN INTERACTIONS 

Let us consider the model Hamiltonian to study the dynamics of S with Hamiltonian Hs, in interaction with Hb- 
Assuming a linear coupling, we have 

Htot — Hs + Hb + Hi 

= ^^^ + H'^^"a«a + X1-9a(4^ + ^^«a), (1) 

A A 

where the operator L acts on the Hilbert space of the system and a\ , a\ are the annihilation and creation operators 
on the environment Hilbert space. The g'^s are the coupling constants that can be taken as real numbers, and the 
uj'^s are the frequencies of the harmonic oscillators that constitute the environment. 

Along the paper we will assume the weak coupling approximation, which consists in considering that the system 
interaction Hamiltonian is g times smaller than the free Hamiltonian, which includes both the Hamiltonian of the 
system and the environment. We will also use for simplicity units in which H — 1. 

The linear model of interaction Hamiltonian is good enough to characterize most of the interactions in quantum 
optics related to processes in which only one photon is involved, and which can be described within the dipolar 
approximation. 

The dynamics of the quantum mean values and MTCF can be computed in two distinct ways: First, by evolving 
vectors on the system Hilbert space (see for instance 13, HI)- In this evolution there are two important quantities 



that describe the action of the environment into the system: the environmental correlation function 

a(i-r)= V5^--^(*-^\ (2) 



A 



which as noted before decays within a time Tc and describes the dissipation of the system, and a certain quantity 
Zt, such that a{t — r) = M[ztz*], where M[] represents a Gaussian average. The function zt can be interpreted as 
a Gaussian colored noise, and describes the stochastic effect of the many environmental degrees on the system wave 
13, 3, [13, [3, II3| ■ Second, in the Heisenberg picture, by evolving system operators. In these equations, the 



vector 



action of the environment is described entirely through the correlation function 
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Note that in the Markov case, a{t — r) « S{t — t) (corresponding to Tc ~ 0), and the quantity zt may represent a 
Gaussian white noise. 

In this paper, we need to compute some two-time correlations of system observables. We wih do it within the 
Heisenberg picture scheme, by using the system of equations derived in l^, 11| in the weak couphng hmit. 



To compute the environment correlation function, it is necessary to perform the sum appearing in ^ , which can in 
principle be done once the value of gx is known for the particular interaction. For instance, in the our case of an atom 
dipolarly coupled with the radiation field, the coupling constants are known, and the quantum number is A = {k, a}, 
where k is the quantum number corresponding to the wave vector, and a corresponds to the two polarization modes. 

However, there are certain problems where the functional form g\ is not known, for instance atoms in contact with 
the phonons in a solid. In those cases, it is often useful to express the sum appearing in ([2|) as 



a{t-T) ^ I dujJ{uj)e-"^^*"-\ (3) 

mm) 



where we have defined the so-called spectral function as J(w) = tt g^d juj — lu\). This spectral function can be 
approximated with some phenomenological models (see for instance 

The spectral function can be expressed as J(w) = G'^ (uj) p{uj) , where G{llj) is the function g\ in the continuum and 
p{uj) is the photonic DOS of the environment. Inserting this relation in ([3]) it is easy to see the important role of the 
photonic DOS in the correlation function, and particularly in its decaying time Tc. 

Finally, we note that despite the weak coupling approximation here used is only valid for weak and medium 
couplings, and therefore it might not capture completely the non-Markovian behavior of the system, it does describe 
some important memory effects in its dynamics. Particularly, the evolution of the system observables depend on their 
past history up to the time Tc that the environment takes to relax back to equilibrium after the interaction. 



ATOMIC EMISSION SPECTRA 

In this section we present the formula necessary to obtain the emission spectra of an atom with non-Markovian 
interaction. To this end, we follow a well known photo-detection model of experiment, the gedanken spectrum analyzer 
that can be found in and which provides an operational definition of the spectral profile. In this scheme, we 

consider that the radiation field emitted by a two-level atom is detected by another two-level atom with frequency 
uj, which is initially prepared in its ground state | 1). The detecting atom has Hamiltonian Hx> = ujaz/2 , and it 
is placed inside a shutter, that only opens during a certain observation time T, during which it receives the emitted 
radiation, being eventually excited. The Hamiltonian of the emitting atom (with levels |1) and |2)) is given by 

Hs = ^('^22 - o-ii) = — 0-z, (4) 

and the total Hamiltonian of emitting atom and radiation field is described by a Hamiltonian Hji, given by Hji = 
Hs + Hb + Y.x9x{L^ax + a{L). 

The Hamiltonian of the total system (detector atom, emitting atom and radiation field) reads as follows, 

H = Hv + Hr + W. (5) 

Here the coupling between the detecting atom Hx> with the radiation field, with Hamiltonian Hr, is dipolar and given 
by a Hamiltonian W , which in interaction image with respect to the detector is given by 

W{t) = Uid^ • E(+) (r, i)e*"* + cri2d^ • E^") (r, t)e-^'^*j , (6) 

where we have considered d^^d^ = df^d-^ = (1 | D | 2) = d^. The super-index T) reminds that these are the 
components of the detector's dipole, and D is the dipolar operator. It is important to note here that the field 
operators E(+) and E'~), correspond to the radiative atoms and the radiation field with Hr. The positive part of the 
field placed in r is defined as 

E^+\r,ro,t) ^J2exAx{r)ax{ro,t)ex (7) 

A 

and E(~)(r,ro,t) = [E'+^ (r, Tq, t)]^ ^2d\. In the last expression (and from now on) we have added explicitly the 



dependency on the position Tq of the source dipole that originates the field. The quantity ex = y with v the 
quantization volume. In terms of the coupling strengths we find that g\ = gxir) — exAx{'r)d ■ ex- 
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As a condition, the time of observation T is much bigger that the inverse of the natural width F of the detecting 
atom excited level. In addition, the Rabi frequency of the emitting atom has to be bigger than the inverse of T. With 
this set up, the idea is to calculate the spectral distribution of the fluorescence light, P{u>,T). This is defined as the 
probability of excitation of the detecting atom at the time of observation T, i.e. 

P{uj,T)^TrR^-D{\2){2\plT)), (8) 

where p{T) is the density matrix of the total system at time T . In the interaction representation, this density matrix 



p{T) = p(0) - I 



dt / dt' 



W{t), W{t'),p{Q) 



(9) 



for an expansion in which p{T) « /o(0). 

Replacing © in ([5]), we get the following expression for P{uj,T), 



P{uj, T) - Trn dt J* dt'e'^^'-^'^d^ ■ E(+) (r, ro, t)pniP ■ E^-) (r, ro, t') 
dt J* dt'e'^^''-'^d^ ■ E(+^(r, ro, t')pRd^ ■ E^-\r, ro, t) j . 



(10) 



If we now change t into t' , and t' into t in the second integral. 



P{cj,T)^TrR |^^^dt^*di'e-(*-*')d^ •E(+)(r,ro,i)p^d^ •E(-)(r,ro,t') 
£dt' dte-(*-*')d^.E(+)(r,ro,tWd^-E(-)(r,ro,0^ , 



(11) 



both integrand becomes the same. On the other hand, the regions of integration are the complementary parts of a 
square of side T, so that the last expression can be written as 



P{u,T) 



dt 



diV-(*-*')g(i)(r,ro;t,i'), 



where the average 



TrR (pji ■ • • ) , and we have defined 

g^'^ (r, ro; t, t') = (d^ • E(-)(r, ro, t)d^ ■ E(+)(r, ro, t')), 



(12) 



(13) 



as the first order correlation of the projection of the emitted field in the direction of the dipole. 

In order to express the last equation in terms of observables of the emitting atoms, it is necessary to express the 
field operator, which is defined in ([7]) in terms of the atomic dipolar momentum. We insert in ^ the solution of the 
Heisenberg equation for aA(ro,<) with the total Hamiltonian ([5]) , i.e. da\{ro,t)/dt = —i[H,a\{rQ,t)], 

a„(ro, t) = aA(ro, 0)6-^'^^^ - i f dTgxL{t')e-'^-^'-^^ (14) 

Jo 

Here, the coupling constant between emitting atoms and surrounding field is gx = gx{vo) — eA^A(ro)d • e^. After this 
replacement, the component of E^+^ in the direction and position of the detector is 



d^.E(+)(r,ro,t) = -*^5Ag? T drL(r) 



(i-r) 



(15) 



where we have defined gf = 5^(r) — e\A\{r)d^ ■ e\ as the coupling constant of the detector with the emitted field. 
In expression we only keep the part of the field corresponding to the radiated field from the dipole (also "source 
field"). The other part, corresponding to the quantum vacuum field 



E[,+'(r,ro,t) = ^eAAA(r)aA(ro,0)e-''^-*e;„ 



(16) 
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does not contribute to photo-detection signals, since provided that the field is in the vacuum state |0), 

d^.E(-'(r,ro,T)) = 



(d^.EW(r,,ro,0)) |0) = (17) 

in p^ . We now introduce a distribution function p{ro), which describes the density of radiating atoms at a given 
positions of the crystal. Assuming that the atomic distribution is the same for each Wigner Seitz cell (WSC) of the 
crystal [25j . we perform an average over the atomic distribution within the crystal, obtaining the following expression 
for the positive part of the radiation field, 

(d^ .E(+)(r,ro,t))ro = ^ / drop(ro)d^ .E(+)(r,ro,i) (18) 

^^e JWSC 

where Nf, = J^rgQ rfro/9(ro) is the number of active atoms in a cell. Considering one single atom at position in each 
cell, such that p{ro) = S{ro — r^) and = 1. Then we have 



(d^ •EW(r,ro,t)}ro = / dr5(r,r„t,r)L(r) 

Jo 



(19) 



with the 

<S(r, r„ t,T)=J2 3?3Ae-^"^(*-^)e'(■^-"■")^ (20) 

A 

where for simplicity in the notation we set gx = g\{ra). Replacing ([19]) in the integrand of ([T2)) . we get 
g(i)(r,r„,i,t') = (g^^) (r, ro, i, t'))ro 

'(d^ . E(-)(r, ra,t))^„) f (d^ . E(+)(r, r„ i'))r„) ) 

(21) 



f dT f dT'S*{v,VaXT)S{r,ra,t',T'){L\T)L{T')) 
JQ Jo 



Following the spectra has the following form. 



p(^,r) = r dt r dt'e^^'-''^g'^^Hv,va.t,t') 

Jo Jo 

i-T pT ( ft ft' 



[ dt [ di'e'"(*-*'4 [ dT [ dT'S*{Y,Ya,t,T)S{v,Va,t',T') 

Jo Jo I Jo Jo 

x{LHr)LiT'))}. (22) 



In order to compute the fluorescence spectra, the limit of T ^ oo has to be taken, so that the signal is observed in 
the stationary limit. In that case, the above formula corresponds to two Laplace transform of a convolution, 

P{u) ^ S*{r,ra,^co)S{r,ra,iu){L\-Lu)L{Lo)). (23) 

Since 

P(^, ^) = P{uj) = Ttr.v (I 2)(2 I p{^)) (24) 
we are dealing with the detector reduced density matrix in the stationary limit. 

Special case: No spatial dependency 

When there is no spatial dependency of the functions 5*(r, ra, i, t), i.e. 

g*k.(ro-r„) _ (25) 
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then 

5*(r,r„t,r) =a(t-T) (26) 
they are the usual correlation functions a{t — t). In that case, we can write (|22p as 

P{w,T) = I dt ( dt'e'^^^'-^'A I dT I dT'a*{t~T)a{t' -t'){L\t)L{t'))\. (27) 



Jo Jo wo ^0 J 

Hence, the fluorescence spectra can be written as 

Piiv) = a*i~uj)aiuj){L\~uj)L{u;)). (28) 

In addition, in the Markov case, such correlation is a delta function, a{t — t) = TS{t — r), and the last formula is 
just 

P(w,T) = r2 r dt r dt'e"^^*-*">{L^t)L{t')), (29) 
Jo Jo 

which in the equilibrium, and with the observation time T —>■ oo, gives rise to the usual expression for the power 
spectra, 

/"OO 

P{uj)=r^ ciTe''""(it(0)i(r)). (30) 
Jo 

The system corr clations {L''{0)L{t)) can be computed with the Quantum Regression Theorem. 

TWO-TIME CORRELATION FUNCTION OF SYSTEM OBSERVABLES IN THE WEAK COUPLING 

LIMIT 

In order to evaluate the spectra from eg. ([^5]) . it is necessary to obtain the evolution of two-time correlations of 
the system observables. A theory to compute multiple time correlation functions in the non-Markovian regime has 
been developed in l^, 11 1. In this section, we briefly present some of the main results, and the formula to compute 







two-time correlations in the weak coupling limit needed in eq. ()23p . 

The evolution equation of the quantum mean value of A, provided that the total initial state is | ^'o) =| V'o) I 0), is 
equal to 

-^(*o I A{h) I *o) = *(*o I [Hsih),A{h)] I *o) 
at I 

+ f ' draiti ~ t)(*o I [LHti), A{ti)]L{T) \ *o) 

dTa*{t, - t){^o I L^{T)[A{ti),L{t,)] \ vPo). (31) 



In the same manner the evolution of the quantum mean value {A(ti)B{t2)) is given by 

d{^o I A{t,)Bit2) I ^o) 

— = zCfo [Hs[ti),A(ti)\B(t2) Wo) 

dti 

rti 

+ / dTa*{ti - t)(*o I L\r)[A{ti), L{ti)]B{t2) \ ^o) 
Jo 

ti 



+ / drait, - t)(*o I [L^{ti),A{ti)]L{T)B{t2) \ *o) 



+ / dTaiti~T){'^o\[LT{ti),A{ti)]B{t2)L{T)\'^o)- (32) 







Equations pip and (|32p represent the evolution of quantum mean values and two-time correlations respectively, 
obtained without the use of any approximation. However, it is clear that these equations are open, in the sense that 



7 



quantum mean values depend on two-time correlations, while two-time correlations depend on three time correlations. 
In general when no approximations are made, TV-time correlation depend on N + 1-time correlations, what gives rise 
to a hierarchy structure of MTCF as described in [l^l ■ 

For practical computations, it is necessary to make some assumptions. In our case we perform a perturbative 
expansion on the parameter g, so that the former equations become closed. Particularly, the perturbative hypothesis 
is applied to the operators L{t) and ^^(t) that appears in equations (PT|) and (|32|. The idea is to transform groups of 
operators with a two-time dependency on groups with a one time-dependency |36[. For instance, the first dissipative 
term in ((3T|) and ([32|) is of type 



L\t) {[L,A]} iU) ^Uf\hO)L{T,U)[L,A]Ui{UO) 
= {Vr-uL^L,A]}iU) + Oig). (33) 

Here, we have used a Liouville super-operator Vt which acts only on the next system operator as VtL = e'^^*Le~'^^'. 
In the same way, the second dissipative term in (|32p can be expanded as 

LiT)B{U+i)=U^\U+iO)LiT,U+i)BUi{U+iO) 
= {Vr^u^,LB}{U+,) + 0{g), (34) 

up to zero order. Note that just as before, the dissipative terms are already of second order due to the appearance of 
the second order function a(t), so that only zero order expansions of the rest of the quantities are needed. Then, up 
to second order in g, we obtain the following equation for quantum mean values, 

-|-(*o I A{ti) I *o) =«(^'o I {[Hs,A]}{ti) I *o) 
at I 

'd™*(ii-r)(*o I {Vr-t,L^[A,L]}{t,)\^o) 
tl 

+ / dTa{ti-T){^o I {[L\A]Vr^uL}{ti) \ V^o) 
Jo 

+0{g'), (35) 
and for two-time correlations 

-^(^0 I A{t,)B{t2) I *o) =*(fo I {[Hs,A]}{t,)B{t2) \ *o) 
at I 

tl 

d™*(ii -t)(*o I {Vr-uL^[A,L]} {h)B{t2) I *o) 

tl 

d™(ii -r)(^'o I {[L\A]Vr-tiL) {h)B{t2) \ V-o) 

dTa(ti - r)(vI/o I {[L\A]} (t,) {[B, Vr^t.L]} (^2) | *o) + O(ff'). (36) 

While the first two terms of ([36)1 are analogous to those of |35)) . the equation for two-time correlations contains 
an additional term that does not appear in the evolution of quantum mean values. Since the Quantum Regression 
Theorem (QRT) states that the evolution of two-time correlations has the same coefficients as the evolution of one- 
time correlations, these equations show that the theorem is not satisfied for non-Markovian interactions. The validity 
of the QRT have been previously discussed by some authors [2I, 23, 23,[24|. Notice that for Markovian interactions 
the last term of (|36[) vanishes, since the corresponding correlation function a{ti — t) = T5{ti — r) is zero in the domain 
of integration from to ^2 , so that the QRT is valid in this case. Note that an equation similar to ([36]) can be derived 
for initial condition {A{ti)B(ti)) , and for t2 > ti- For the Markov case, it can also be shown that this equation is 
equivalent to the one given by the QRT for such an initial condition. For our purposes, we shall use equation to 
evaluate the evolution of non-Markovian two-time correlations. 



FLUORESCENCE SPECTRA OF A TWO LEVEL ATOM 



The formulas derived up to now are valid for the emission spectra of an arbitrary atom, provided that the interaction 
Hamiltonian has the linear form ([T]) and that the weak coupling approximation is valid. In this section, we show how 
to express the Hamiltonian for a driven two-level atom, so that a particular form for the coupling operator L is given. 



8 



As in [26[, we consider the interaction Hamiltonian Hsl of the atom with the classical coherent monochromatic 
laser field in the usual rotating wave approximation form 

Hsl = e(a2ie-*('"^*+'^^) + (Ti2e^("^*+"*^)), (37) 

where lol is the frequency of the laser in the coherent state ae"*"^*, with a — jaje""^^. The quantity e = d2ie is 
the Rabi frequency, and |e| = lo / 2eQV \/ Nl^l the laser field magnitude, where v is the volume of the cavity, and 
Nl =\ \ a-iid cl, are respectively the mean number of photons and the polarization of the laser mode. The phase 
(f)T is a global phase defined as ipT = <t>L — 7r/2. Because of the magnitude of the laser field, the Hamiltonian Hsl 
should be considered as part of the non- interacting Hamiltonian Hq, where in this case Hq = Hs + Hb + Hsl- We 
can eliminate the explicit dependence on the laser frequency by changing to a rotating frame with a frequency o^l by 
using the unitary operator 

If Aa ^ uj\ — lul and A^^ — lus — ojl the rotated Hamiltonian H' = Hq + iJj, where 

HQ = ^Axalax + l-uJs(J3 + e[(j2i+(Ti2] (39) 



2 

A 



and 



H'j = i'^gx{ai2a{- ax<72i), (40) 

A 

can still be expressed in a simpler way by projecting it into the dressed atomic basis. The new Hamiltonian H = 
V-^H'V, where 

will be of the form ([1]). The constants appearing in the unitary transformation matrix V are c — cos0 and s = 
sint/), where the angle, is given by sin^ </> = ^[1 — sff'^(^SL)/\/c/A|^J + l]. The non interacting dressed state 
Hamiltonian Hq — Hs + Hb + Hsl is equal to, 

Ho^nR3 + Y,^xa{ax, (42) 

A 

and the interaction Hamiltonian Hj has the same form as in ([T]) once the interaction operator L is defined as: 

L = csR3 + c^Ri2-s^R2i. (43) 

Here, Rij — are the atomic operators defined in the dressed state basis {|1), |2)}, and i?3 = R22 — Rii- The 

quantity fl — [e^ + A|^/4]-'^/^ is called the generalized Rabi frequency. 

According to ((4T|) . the relation between the dressed atomic operators and the bare atomic operators is the following, 

(T12 = CSi?3 + C^Ri2 - s2i?21 
(721 = CSi?3 - S^i?i2 + C^i?21 
(T3 = (c2-s2)i?3-2cs(i?i2+i?2l). (44) 

It can be easily verified that the case of spontaneous emission corresponds to switching off the laser, and therefore 

L = (712. 

Once the coupling operator L is known, the two time correlation {L^ {t)L{t')) can be easily computed with equation 
([36|) . Now, in order to obtain the emission spectra with (|23|) . it is necessary to derive the Laplace transform of the 
spatial dependent correlation function, 5* (r, r^, w). Notice that the interaction with the field enters into the formula 
(j23p in two ways: 



Indirectly, through the correlation function a{ti — r) entering in equation ([361 
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• Directly, through the function S* {r,ra,uj). 

Let us make here a remark. In order to compute the atomic dynamical equations, we assume that the field is composed 
of extended modes, in the sense that they correspond to a real wave vector. Hence, the spatial dependency is assumed 
to be exp (ik • r) w 1 within the optical region. This assumption that the field is composed only of extended modes 
corresponds to consider the photonic crystal as an infinite and perfect structure, so that evanescent modes are in 
principle not present. In other words, even if evanescent modes are solutions of the eigenvalue problem, normally they 
cannot be excited inside a perfect (infinite and without defects) photonic crystal. This is because the do not satisfy 
the translational simmetry boundary condition of the crystal. 

Nevertheless, a defect placed in the crystal (in the our case an emitting atom) may substain such modes. To be 
more specific, once an impurity atom with resonant frequency within the gap region is placed within the crystal, 
its emission will be in the form of evanescent modes (i.e. with imaginary wave vector). Then, the approximation 
exp (zk • r) 1 is no longer valid for such field. As will be explained in detail later, this fact is very relevant to 
compute the emission spectra, which appears to be strongly dependent on the atom-detector distance through the 
spatially dependent function 5*(r,ra,w). 

In the next section, we present the calculus of a{t — T) and 5*(r, r^, uj) necessary to compute the system fluctuations 
and the fluorescence spectra respectively. Both quantities depend on the particular interaction considered. We study 
them for the case of an atom interacting with the radiation field within a photonic crystal. 



ENVIRONMENT CORRELATION FUNCTIONS FOR AN ATOM IN A PEG 

In this section, we calculate the correlation function characteristic of the dipolar coupling of a two level atom to 
the modified radiation field within a PC. In that case, the coupling constants entering in ([2]) are 



5A = -^^^-i2e...d,,e--, (45) 

where di2 is the atom dipolar moment, Bko- is the unit vector in the direction of the polarization a for a given wave 
vector k, and eo is the electric vacuum permittivity. In the last expression we have already expanded the potential 
vector in terms of plane waves. This is also valid for radiation fields in photonic crystals, but with the additional 
requirement that these plane waves are periodic, i.e. satisfy the Bloch theorem. As we have already noted in the 
introduction, the presence of a periodic dielectric structure produces the appearance of optical bands, which are 
themselves periodic structures in the reciprocal lattice. Periodic eigenvalues (the frequencies), correspond to periodic 
eigenvectors (the field). However, although the quantity A = {k, cr, n}, where n denotes each band, this index will be 
here dropped, since we consider the dynamics of the two level system with rotating frequency nearby (or inside) a 
single band. 

In order to calculate a{t), and S{r,ra,t), the sums appearing in equations ^ and ([20ll respectively, have to be 
performed, by using the functions g\ as defined in (|45p . In that way we have 

a{r) = 7(f )^ E / dk^^i:^e--e^)^ (46) 

„ JlBZ '^W 



and 

where is the unitary vector corresponding to the dipolar moment of the detector, and the integrals are performed 
over the first Brillouin zones of the crystal. Also, it is clear from that the constant 7 — (jj^2'^i2/(2£ofi-)- In order 
to calculate the integrals appearing in (pS)) and (|T7|) . it is clear that the dispersion relation w(k) is needed. 

The dynamics near the edge of the band are often described through an effective mass a ppr oximation to the full 
dispersion relation, based on its expansion in the vicinity of the band edge (see for instance |26l. [27. 28 . 2^). Within 



this approximation, and considering the simplest case of a cubic lattice, the dispersion relation has the parabolic 
form w(k) = ujc + A{\fi — ko)^, where ko is the origin of the first Brillouin zone of the crystal (which is the unitary 
cell in k space) about which we perform the expansion in each direction, Uc is the frequency of the band edge, and 
^ is a constant that depends on the specific photonic crystal considered. Although it is valid for the description 
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of several important physical phenomena occurring in PBG materials 0, H^, [s^, 31 , [s^l , this model of dispersion 
relation becomes inaccurate for a good description of the short time dynamics. The dispersion relation proposed in 
[ist overcomes this problem by reproducing the periodic structure of the bands, 

u}{k) = A + — {cofi{kxa) + cos{kya) + cos{kza)). (48) 
o 

Notice that this function reduces to the parabolic model for k « kg. For a three dimensional band structure with 
the periodic dispersion relation (|48p . the correlation function can be obtained from the equation (|46p once it has 
been assumed that the function \ek.cr ■ Wrfp/a;(k) is a slowly varying function in the IBZ. The resulting integral is 
analytical, and gives rise to the following result: 



aMr) = e J,{ — ) 

= gloe^'^'jU^), (49) 

where Jq is the zero order Bessel function. 

Let us now consider the spatial dependent correlation function, iS(r, r^, t, t). As before, we assume that the quantity 

^'^'' '^ "^^{k) '" ^ ^ slowly varying function in the reciprocal space (particularly in the region nearby the symmetric 
point ko), in therefore we can express the function (|47p as 

5(r,r„,T) ^7( — )3y l^ko,^ ' "rf||gko,. ■ tl^l f ^^^^^u.(l.)r ^^ir-r^)■l. ^ (gp) 

^ w(ko) JiBZ 

Its Laplace transform is then given by 



27r ^ ko JiBZ «wk)-w 



In order to perform the last integrals analytically, it is convenient to consider the parabolic dispersion relation and 
the limit in which the detector and the emitting atom are far away from each other. 

In that way, we can make a change of variable q = k — kg, so that we have w(k) — uJc + Aq^ and extend the limits 
of integration to infinity. Then 

S{r,ra,Lo) ^ -^Q d^q^. ^ -Q-J / dqq^- (52) 

J J J UJ[q)-uj id Jo uj{q)~uj 

We have also assumed that q = q(sin6'cos^, sin6'sin0,cos6'), and the vector d = r — being parallel to the z axis. 
The constant 

Q - 7(|^)^e--'^(l - - ^^^^)' (53) 

provided that the c?]^ — di2- Considering that 9 and 9^ arc the angles between ko and and respectively, we 
have simply 

Q = ^{-^fe'^o'^sin^eBin^e'^ (54) 

Notice that in our case, since we are dealing with a cubic lattice, only one symmetric point ko and the dispersion 
relation around it is being considered. In a more general case, the band edge is associated with a finite collection of 
symmetry-related point kg, so that the calculus has to be taken into account the dispersion relation near each of the 
points, w(k) — uoc + A{k — kp)^. In such a case, the constant 

Q^^i^-^fYe^o'^siu'e.sinHf (55) 



Performing the last integral of ([5^ . we just obtain 



5(r,r„c.) = ^e-'*/' (56) 
id A 
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where I 

As displayed in Figure H]), an increasing distance in frequencies with respect to the band edge (placed at Wc = 0), 
gives rise to a decreasing localization length I, and therefore to a decreasing value of the function S [sF 




FIG. 1: Real and Imaginary parts of the function (|56|) . where we have chosen ljc = 0, and A — 1.. The localization length is 
also displayed for the same parameters, showing a decreasing for increasing detunings from the band edge. 



Numerical results 



In this section we compute the florescence spectra of a two-level atom placed as an impurity in a photonic crystal. 
To this order, we will study two cases, the florescence spectra when no spatial dependency is considered, equation 
and when the spatial dependency is considered (i.e. the distance from the emitting atom to the detector), in 
equation ((23)) . In both cases, we need to compute the Laplace transform of the system fluctuations (L^(T)L(r')), 
obtained with equation ([36]) by using the coupling operator L given by ([43]) . We show that the spatial dependency 
should be considered when part of the radiation is emitted in the form of evanescent waves, since in that case the 
approximation exp (ik • r) w 1 is no longer appropriate. Evanescent modes are emitted for frequencies within the gap. 



No spatial dependency considered 



In this section we compute the atomic fluorescence spectra according to equation (I28p . by using the correlation 
function As we will see in the next section, the formula ([25]) is suitable to compute the emission spectra 

corresponding to frequencies outside of the gap. 

Figure ^ shows a case in which we have chosen the Rabi and the laser frequencies {ft and tj^) in such a way that 
the we have placed the three components of the MoUow triplet inside the band. This corresponds to the setup of 
Figure ([2]-a)). The appearance of a three component fluorescence spectra for a driven atom is due to the dressing of 
the atomic levels by the laser Q . 

In Figure (|4]) we show a case in which the two side-bands of the triplet are placed inside the gap (See Fig. (HJ-b)). 
As we will clearly show in the next section, this two lines corresponds to emission of evanescent modes, which have 
purely imaginary wave vectors. In consequence, the approximation exp (zk • r) « 1 over estimates their magnitude, 
and a more accurate description of the spatial dependency is needed. Only with such description we will be able to 
account for the fact that evanescent modes should not be detected in the far field. This is done in the next section. 
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FIG. 2: We display the difFerent arrangements we consider. Figures a) and b) corresponds to the two setup described in Section 
0, while Fig. c) corresponds to the setup of Section (). 




FIG. 3: A comparison between the MoUow triplet obtained with Markovian interaction and the triplet obtained with the 
non-Markovian interaction within a PC is here shown. The tree components of the triplet are placed within the band, that 
corresponds to frequencies A — B < uj < A + B, where A — 1 and B — 1. The assymetry between the left and the right 
side-bands described in [s^l is also described with our model. 



Spatial dependency considered 

From the last section we have concluded that for a correct description of the emission from the gap, it is very 
important to take into account the spatial dependency, i.e. the distance of the detector from the emitting atom. To 
this order, we consider the equation (|23p. and the Laplace transform of the spatial correlation function (|56p that is 
obtained when considering the parabolic dispersion relation a;(k) = u;c+^(k— ko)^. Moreover, in order to compute the 
system fluctuations (f36l) . we take the (not spatially dependent) correlation function corresponding to this dispersion 
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FIG. 4: In contrast with Fig. ([Sjl, now we place the left and right emission lines within the gap. When the distance between 
atom and detector is not considered, there is no filtering of the fluorescent emission that occurs within the gap. 



relation, which according to [26[ is 



This correlation is just ((49|) taking the long time limit, so that the constants can be related as [l3| 

r,l/2 ^ ffio /^n3/2 

8 ^b' ' 



(57) 



(58) 



When using the parabolic dispersion relation, only two frequency regions have to be taken into account: a gap, 
for oj < u!c, and a band for u> > ujc- We chose our parameters {Aal — 0, ^^l — 1, and il) such that the left MoUow 
component is placed in the gap, while central and right component have frequencies within the band (See Fig. (QJ-c) 
for a schematic view.). In figure ([5]) we show the fluorescence spectra for a certain distance d of the detector (expressed 
by taking the periodicity of the crystal as the unity) which is small enough as to be still detecting some evanescent 
modes from the gap. Different choices of fl shows that the more deep in the gap the emission occurs, the less we 
are efficient to detect it, due to the decreasing of the localization length. Figure (O shows the spectra for a fixed 
value of but varying the distance of the detector, in such a way that for d large enough, a complete suppression of 
fluorescence from the gap (corresponding to the left side-band) occurs. 

Concerning this, two important remarks have to be made here. 



In contrast to 3j| and [33|, where a filter Q{uj — uJc) is introduced for emission frequencies within the gap, in 
the description presented here no filtering effect occurs unless the spatial dependency is considered, and the 
detector is placed in the far field. 

A second remark concerns the difference between the spontaneous and the fluorescence spectra. Even considering 
the spatial dependency and the detector in the far field, some radiation is detected within the gap region in the 
spontaneous emission spectra. The reason is that there is a transient field, called diffusive field, that appears 
at the beginning of the atomic evolution. Since this field can be emitted in the gap frequencies, it carries away 
energy from the atom. Hence, it is responsible of the slight dissipation that occurs in the atomic population 
of the excited level at the initial time steps of the evolution [l5|, |26|. A description of this field, including its 
spatial dependency is found in [sH]. In this paper, this is done by solving the atom dynamics within the one 
photon sector [3§]. 

In contrast, the fluorescence spectra deals with the signal emitted by the atom in the stationary limit, and 
therefore the transient held is no longer detected. For that reason, the radiation within the gap is only composed 
of evanescent modes, and cannot be observed when placing the detector far away from the atom. 
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FIG. 5: Power spectra d?P{uj) for different values of the Rabi frequency f2, and a fixed distance of tlie detector d = 10. In all 
the cases the left sideband of the MoUow triplet is emitted in gap frequencies. While the distance d is small enough to detect 
the evanescent modes, the fact that the localization length is smaller for increasing distance to the band edge uic ~ 0, makes 
the left MoUow component be increasingly small for increasing values of Q. 

In summary, it is true that the fihcring appears as a consequence of the special form of the dispersion relation 
of the field in the PC. However, we show here that this filtering only occurs provided that the signal is in the 
stationary limit, and the distance between the emitting source and the detector is large. 

CONCLUSIONS 

We present in this paper a general formula to compute the emission spectra of an atom linearly coupled to an 
environment with non-Markovian interaction, and within the weak coupling approximation. This formula is explicitly 
dependent on the distance between the emitting atom and the detector, d. We apply it to compute the fluorescence 
spectra of a two-level atom in a photonic crystal. In those materials, the electromagnetic field have very specific 
features, showing a gap or range of frequencies where the propagating modes are not allowed. In that frequency region, 
the emission spectra depends critically on d, since the emission occurs in the form of evanescent (non propagating) 
modes, which decay exponentially from the atom. For that reason, we show that in the far field these modes are not 
detected. We note that the description of this evanescent field can only be made with a non-Markovian description, 
since a Markov theory, with dissipation rates calculated through the Fermy Golden Rule, predicts no emission at all 
within the gap region. 
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FIG. 6: Power spectra cPP{lo) for different values of the the detector distance d, and fixed Rabi frequency = 0.55, so that 
the left sideband of the MoUow triplet is emitted in gap frequencies. It is observed that for distances large enough this line is 
no longer detected, an a complete filtering of the spectra within the gap is obtained. 
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